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Abstract. We used 1627 faint (15.5 < -R < 23) stars in five fields of the Calar Alto Deep Imaging Survey (CADIS) to estimate 
the structure parameters of the Galaxy. The results were derived by applying two complementary methods: first by fitting the 
density distribution function to the measured density of stars perpendicular to the Galactic plane, and second by modelling the 
observed colors and apparent magnitudes of the stars in the field, using Monte Carlo simulations. The best-fitting model of the 
Galaxy is then determined by minimising the C-statistic, a modified x^- Our model includes a double exponential for the stellar 
disk with scaleheights hi and h2 and a power law halo with exponent a. 24 480 different parameter combinations have been 
simulated. Both methods yield consistent results: the best fitting parameter combination is a = 3.0 (or a = 2.5, if we allow 
for a flattening of the halo with an axial ratio of (c/a) — 0.6), hi = 300 pc, /i2 ~ 900 pc, and the contribution of thick disk 
stars to the disk stars in the solar neighbourhood is found to be between 4 and 10 %. 
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1. Introduction 

The stellar structure of our own Galaxy has been studied in- 
tensively for almost 400 years now, and although, or perhaps 
exactly because of the fact that the subject to be analysed is 
our own neighbourhood and surrounds us, there are still a large 
number of unknowns and the exact structure parameters are 
still a topic of debate. 

In the standard model dBahcall & Soneiral 1 1980a b) the 
Galaxy is of Hubble type Sbc, consisting of an exponential 
disk, a central bulge, and a spherical halo. The vertical struc- 
ture of the disk follows an exponentia l law (p oc 6"^^^ "= ) with 
scaleheight hz — 325 pc. However, as lGiknore&R'eidl(ll983h 
showed, the data can be fitted much better by a superposition 
of two expone ntials with sca leheights hi of 90 to 325 pc and 
/i2 ~ 1300 pc ( lGilmorell984 t). It is not clear whether this devi- 
ation from the single exponential is due to a distinct population 
of stars, although this is s uggested by the different kinemat- 
ics and lower metallicities ( lFreemanlll992l) . it is referred to as 
"thick disk" in the literature. 

In the last two decades highly efficient, large area sur- 
veys have been carried out, and large quantities of high-quality 
imaging data have been produced. The mere existance of a 
thick disk component is fairly established today, however, the 
size of its scaleheights are still under discussion ('N orriJl999t 
[sieg el et al. 2002), as several authors claim to have found a 
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considerably smaller value for the scaleheig ht of the thick disk 
than the canonical one of h2 ~ 1.3 kpc: iRobin et alJ tOOdi 
found 750 pc, and lChen etatl ( 1200 ll) found the scaleheight to 
be between 580 and 750 pc. 

In general there are two different approaches to deduce 
Galactic Structure: the B aconian, and the Cartesian ansatz, as 
iGilmore & Wvs3 lll987l) call it. 

The Baconian ansatz tries to manage with a minimum of 
assumptions: by a pure measurement of the stellar distribution 
function by means of starcounts. Distances are estimated for 
each star and then their number density in dependence of dis- 
tance from the Galacti c plan e is calculated (G ilmore & Reicjj 
Il983t iReid et aljfl993 Il997t iGould et alJil998l) . The param- 
eters can then be determined by fitting distribution functions 
to the data. In a first paper on CADIS deep star counts 
dPhleps et alJ2000t Paper I in the following) we presented first 
results based on k, 300 faint stars (15.5 < R < 23) in two 
CADIS fields covering an area of ~ 1/15 D° in total. From 
these data we deduced the density distribution of the stars up 
to a distance of about 20 kpc above the Galactic plane, us- 
ing the Baconian ansatz. We found hi — 280 pc, and unam- 
bigously the contribution of the thick disk with /i2 on the order 
of 1300 pc. 

On the other hand, as the main structural features of 
the Milky Way have been identified, it has been possible 
to design synthetic models of the Galactic stellar popula- 
tions (.Bahcall & Soneira.. 1981: .Gilmore. 1984^ .Robin & Crez3 
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1986t iReid & Maiewskil Il993t iMendez & van Altenal fl998l 
Chen et al 2001), assuming the forms of the density dis- 

tribution function f or the different compone nts of the Galaxy. 
This method is what lOiknore & Wvsel ( ll987l) call the Cartesian 
ansatz- 

The exact structure parameters can then be deduced by 
comparing model and star count data. 

In this paper we use both approaches using the stellar com- 
ponent of the CADIS data. Deep {R < 23) multi color data is 
now available for 1627 faint stars in five fields at high Galactic 
latitude and different Galactic longitudes. 

In Section|2]the Calar Alto Deep Imaging Survey is briefly 
described. In Section0]the results we deduced by applying the 
"classic", Baconian method are shown, whereas in Section |5] 
we show the results given by the modelling approach. A sum- 
mary and discussion is given in Section|6] 

2. The Calar Alto Deep Imaging Survey 

The Calar Alto Deep Imaging Survey was established in 1996 
as the extragalactic key project of the Max-Planck Institut fiir 
Astronomic. It combines a very deep emission line survey car- 
ried out with an imaging Fabry-Perot interferometer with a 
deep multicolour survey using three broad-band optical to NIR 
filters and up to thirteen medium-band filters. The combination 
of different observing strategies facilitates not only the detec- 
tion of emission line objects but also the derivation of photo- 
metric spectra of all objects in the fields without performing 
time consuming slit spectroscopy. 

The seven CADIS fields measure w 1/30 0° each (11" x 
11") and are located at high Galactic latitudes to avoid dust 
absorption and reddening. In all fields the total flux on the 
IRAS 100 fim maps is less than 2 MJy/sr which coiTesponds to 
Eb-v < 0.07. Therefore we do not have to apply any colour 
corrections. As a second selection criterium the fields should 
not contain any star brighter than « iQ"^°-9 tjjg CADIS R 
band. In fact the brightest star in the five fields under consider- 
ation has an R magnitude of 15.42™°^. 

All observations were performed on Calar Alto, Spain. In 
the optical wavelength region the focal reducers CAFOS (Calar 
Alto Faint Object Spectrograph) at the 2.2 m telescope and 
MOSCA (Multi Object Spectrograph for Calar Alto) at the 3.5 
m telescope were used. The NIR observations have been car- 
ried out using the Omega Prime camera at the 3.5m telescope. 

In each filter, a set of 5 to 15 individual exposures was 
taken. The images of one set were then bias subtracted, flat- 
fielded and corrected for cosmic ray hits, and then coadded 
to one deep sumframe. This basic data reduction steps were 
done with the MIDAS software package in combination with 
the data reduction and photometry package MPIAPHOT (de- 
veloped by H.-J. Roser and K. Meisenheimer). 

2. 1. Object detection and classification 

Object search is done on the sumframe o f each filter using 
the s ource extractor software SExtractor jBertin & ArnoutsI 
ll996l) . The filter-specific object lists are then merged into 
a master list containing all objects exceeding a minimum 



S/N ratio in any of the bands. Photometry is done us- 
ing the program Evaluate, w hich has been develot^ed by 
Meis enheimer and Roser (see iMeisenheimer & Rosen Il986l 
and Roser & Meisenheimeilll99ll) . Variations in seeing be- 
tween individual exposures are taken into account, in order to 
get accurate colours. Because the photometry is performed on 
individual frames rather than sumframes, realistic estimates of 
the photometric errors can be derived straightforwardly. 

The measured counts are translated into physical fluxes out- 
side the terrestrial atmosphere by using a set of "tertiary" spec- 
trophotometric standard stars which were established in the 
CADIS fields, and which are calibrated with secondary stan- 
dard stars cOke.. 1990; . Walslt 1 995.) in photometric nights. 

From the physical fluxes, magnitudes and colour indices 
(an object's brightness ratio in any two filters, usually given 
in units of magnitudes) can be calcul ated. The CADIS m ag- 
nitude system is d escribed in detail in IWolf et al.lll200Ibll and 
iFried et alJ(l200ll) . 

The CADIS color is defined by: 

(fe-r) =2.51og|| . (1) 

Here is the flux outside the atmosphere in units of Photons 
nm^i in the CADIS filter k. 
b — r can be calibrated to the Johnson-Cousins system (the 
CADIS Rc is very close to the Cousins R) by using Vega as a 
zero point: 

^Vcga(A = 440nm) 



{B-R)c = (6-r) + 2.51og 
= (&-r)+ 0.725 



^{^c.a(A = 648nm) 



(2) 



With a typical seeing of 1'.'5 a morphological star-galaxy 
separation becomes unreliable at i? ^ 21 where already many 
galaxies appear compact. Quasars have point-like appearance, 
and thus can not be distinguished from stars by morphology. 
Therefore a classification scheme was developed, which is 
based solely on template spectral energy distributions (SEDs) 
(Wolf et al. 2001b a). The classification algorithm basically 
compares the observed colours of each object with a colour 
library of known objects. This colour library is assembled from 
observed spectra by synthetic photometry performed on an ac- 
curate representation of the instrumental characteristics used 
by CADIS. As input for the stellar library we used the cata- 
logue by Pickles (1998). 

Using the m inimum variance estimator (for details see 
IWolf etal.l20"01 b). each object is assigned a type (star - QSO - 
galaxy), a redshift (if it is not classified as star), and an SED. 

Note that we do not apply any morphological star/galaxy 
separation or use other criteria. The classification is purely 
spectrophotometric . 

Five CADIS fields have been fully analysed so far (for co- 
ordinates see Table0. We identified 1627 stars with R < 23. 
The number of stars per field is given in TableQ] 

3. The stellar density distribution 

Since none of our fields points towards the Galactic center and 
all of them are located at high Galactic latitudes, neither a bulge 
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Table 1. The Galactic coordinates of the five fields investigated 
so far, and the number of stars per field, R < 23. 



CADIS field 


I 


b 


stars 


Ih 


150° 


-59° 


200 


9h 


175° 


45° 


300 


13h 


335° 


60° 


517 


16h 


85° 


45° 


390 


23 h 


90° 


-43° 


220 



component nor spiral arms are needed to describe the observa- 
tions. Our model of the density distribution of the stars includes 
only a disk and a halo. 

3. 1. The distribution of stars in tine Galactic disk 

The vertical density distribution of stars in the disk is described 
by a sum of two exponentials (the "thin" and a "thick" disk 
component): 



n2 exp - — 
h2 



(3) 



where z is the vertical distance from the Galactic plane, ni 
and n2 are the normalisations at z = 0, and hi and h2 are the 
scaleheights of thin and thick disk, respectively. 

The radial change of the density can be taken into ac- 
count by assuming an exponential decrease with scalelength 
hr- The value of h^ is still poorly determi ned, measurements 
range from values as small as hr = 2 kpc ('Robin etalJll992t 
[Ruphvetal. 1996; Oiha e t ap[l996t Bienavme 19 9^to K = 
3.5 kp c ( iBahcall & Soneiral Il980bril98li IWainscoat et alJ 



For each field with Galactic longitude I and latitude b, re- 
spectively, we can describe the distribution of the stars in the 
disk by 



Pdir,l,b) 



ni exp 



r sin b 
hi 



n2 exp 



• exp 



^ tan h 



COS I + sin b 



rsmb 
h2 

Rq 



The distance of the sun from the Galacic center is assumed to 
be Rq — 8 kpc. 

3.2. The distribution of stars in the Galactic halo 

The space density correspon ding to the r^J ^ law in projection 
has no simple analytic form ('Youna'1976!). There are different 
possibilities to find an appro ximation, the simplest of which is 
a power law dFuchs & JahreiBi 19981) of the form 



R% 



y/2 



with an arbitrarily chosen core radius of Tc 



1 kpc, and 



r^R: 



- 2R, 



tan 6 



cos / + 



(c/ay 



(5) 



(6) 



Several authors found that the stellar halo is significa ntly 
flatte ned, with an axial r atio of (c/a) around 0.6 ( Robin et alJ 
E00(llT,emon etalJ2004 



4. Star Counts - a direct measurement 

Now, as we have more than twice the number of fields available 
than in our first analysis, and due to the increased number of 
filters used for the classification of the objects, the statistics has 
increased considerably, from ^ 300 stars in Paper I to ~ 1630 
now. 

Figure Q shows the distribution of measured colors of the 
selected stars in the field. As in our first analysis, we use this 
bimodal distribution of colors, which comes from the interac- 
tion between luminosity function, density distribution and ap- 
parent magnitude limits (15.5 < R < 23), to discriminate be- 
tween disk and halo stars. The nearby disk stars are intrinsi- 
cally faint and therefore red, whereas the halo stars have to be 
bright enough to be seen at large distances, and are therefore 
blue. According to the different directions of the fields the cut 
is located at slightly different colors, as indicated in Figure ^ 
Unfortunately it is not possible to dicriminate between thin and 
thick disk stars in the same way, so in the calculation of abso- 
lute magnitudes and thus distances we have to treat them in the 
same way. 

We repeated the calculation of the photometric parallaxes 
of the stars exactly as described in Paper I. Absolute mag- 
nitudes are inferred from a mean main sequence relation in 
a color-magnitude diagram (lT,andll992l) . The main sequence 
converted into a vs {b — r) relation (see Paper I) can 
be approximated by a fourth order polynomial in the range 
-2 < (6-r) < 1.8: 



Mr = co + ci{b-r) + 02 {b - rf 
+ C3 {b-rf +Ci {b-rf, 



(7) 



(4) the parameters of which are: 



Co = 4.01236 ci = 4.12575 cz = -1.89076 
C3 = 0.762053 C4 = 0.341384 . 

This main sequence approximation is valid for all stars in 
our sample since we can be sure that it is free from contami- 
nation of any non-main sequence stars, as the faint magnitude 
intervall we observe (16 < R < 23) does not allow the detec- 
tion of a giant star 

The mean main sequence relation is valid strictly only for 
stars with solar metallicities, whereas our sample may contain 
stars spread over a wide range of different metalicities. The in- 
fluence of the varying metalicities on the main sequence color- 
magnitude relation is taken into account as described in Paper 
I: Halo stars are supposed to be metal poor, and are known 
to be fainter than disk stars with the same colors, so the rela- 
tion is shifted towards fainter magnitudes by AMr = 0.75. 
This value is the mean deviation from the mean main sequence 
defined by the CNS4 stars UahreiB & Wielea.1997.) of a sub- 
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sample of 10 halo stars for which absolute R magnitudes were 
available (JahreiB, priv. com).' 

The spread in a two color diagram {b — r) versus (r — i) 
(that is the CADIS color between Rc and /sis, analog to Eq. 
[3 becomes significant at {b — r) « 1.0, see Fig.|2] The maxi- 
mal photometric error for the very faint stars is 0.15™°^. Here 
metallicity effects will distort the relation between the mea- 
sured {b — r) colors and the spectral type (temperature) and 
thus lead to wrong absolute magnitudes, so we have to correct 
for metallicity in order to avoid errors in the photometric par- 
allaxes. 

The Rc filter is strongly affected by metallicity effects like 
absorption bands of Ti02 and VO molecules in the stars' atmo- 
sphere, whereas the Bq and the medium-band filter /gi5 (the 
wavelength of which was chosen in order to avoid absorption 
bands in cool stars) are not. So in a first approximation we can 
assume the "isophotes" of varying metallicity in a (6 — r) ver- 
sus (?■ — i) two color diagram to be straight lines with a slope 
of —1, along of which we project the measured colors with 
{b — r) > 1.0 onto the mean main sequence track which in the 
interval 1.0 < {b — r) < 1.8 is defined by 



(r-t) = 0.39 {b - r)t,,, - 0.36 (6 - r) 
+ 0.09 {b-r) + 0.06 . 



3 

corr 



(8) 



This projection implies that stars with (6 — r)corr 1-8 
cannot exist in Fig. |3l which shows the spatial distribution of 
metallicity corrected (6 — 7')corr colors (the limit is indicated 
by the dashed-dotted line). Both the upper and lower magni- 
tude limits lead to selection effects which have to be taken into 
account. 

4.0.1 . Completeness Correction 

We correct for incompleteness in the following way (for a de- 
tailed derivation of the completeness correction and the cor- 
responding errors we refer the reader to Appendix 1 of Paper 
I): First, we divide the distance in logarithmic bins of 0.2 and 
count the stars up to distance dependent color (that is, lumi- 
nosity) limit, up to which stars can be detected, see Fig. |3] 
The nearest bins (—0.8 < logr < —0.2) are assumed to be 
complete. For the incomplete bins we multiply iteratively with 
a factor given by the ratio of complete to incomplete number 
counts in the previous bin, where the limit for the uncorrected 
counts is defined by the bin currently under examination (j): 



N. 



(9) 



where Nj is the number of stars in bin j, N is the number of 
stars in the previous bin (j — 1), up to the limit given by the bin 
j, and N is the number of stars from that limit up to the limit 
given by bin j — I- 



' Note that this artificial separation may well lead to wrong absolute 
magnitudes for individual stars (e.g. a disk star with r < 1 kpc but 
6 — r < 0.7), but should be correct on average. 



With the poissonian errors (Jn — y/N, CTjv' = and 



(T»f" = V N" the error of the corrected number counts be- 



'N 

comes: 



N. 




(10) 



Note that the completeness correction is done for each field 
separately, thus the different color cuts by which we separate 
disk from halo stars can be taken into account. Also since the 
color-magnitude relation is different for disk and halo stars, we 
calculate the correction factors separately for the two popula- 
tions. The correction factors are usually in the range 1 < C < 
10, but can be as large as 20 for the last bin of the halo. 

With the corrected number counts the density in the loga- 
rithmic spaced volume bins (Vj = ^u}{rj^i — r|)) can then be 
calculated according to 



Pj = 



Tycorr 



(11) 



For every logarithmic distance bin we use the mean height z 
above the Galactic plane < Zj >= sin&- < r >, where < 

logr >= logrj + (logrj+i - logr,:)/2 = logr^ + 0.1 ^< 
r >= 1.259n. 

From the estimated density distribution in the five fields we 
now deduce the structure parameters by fitting the distribution 
function to the data. 



4. 1. The distribution of stars in ttie Galactic disk 

We first calculate the density of the stars which we assume to 
be disk stars, that is, with colors redder than the color cut (see 
Figure[2in the field under consideration. 

As the nearest stars in our fields have still distances of about 
200 pc the normalization at z = has to be established by 
other means: We take stars from the CNS 4 (Fourth Catalogue 
of Nearby Stars, IJahreiB & WielenI ( Il997l) l which are located 
in a sphere with radius 20 pc around the sun. The stars in our 
normalization sample are selected from the CNS4 by their ab- 
solute visual magnitudes, according to the distribution of abso- 
lute magnitudes of the CADIS disk stars (6.5 < My < 14.5). 

Since the thin disk completely dominates the first kpc and 
there is no indication of a c hange of slope within this dis- 
tance ('Fu chs&Wieienlll993l) . it is justified to fit the density 
distribution in this range with a single exponential compo- 
nent (taking of course the radial decrease into account). The 
sum of thick and thin disk stars in the solar neighbourhood, 
{ni + n2) — 6.566 • lO^kpc^"^ with an uncertainty of 10%, 
is given by the CNS4 stars. We estimate the scaleheight hi 
from this fit, assuming two extreme values of the scalelength: 
hr = 2.0 kpc, and hr — 3.5 kpc, respectively (see Table|2j. For 
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-10 12 

log r / kpc 



Fig. 3. Spatial distribution of metallicity corrected (6 — r)coir- The solid line represents the distance dependent upper color limit 
at i? = 23™°^, the dashed line is the lower limit due to the selection criterium of the fields (no star brighter than R w 15.5). Since 
the metal-poor halo stars are intrinsically fainter, the color limits are shifted accordingly. The horizontal line at {b — r)cori = 0.75 
denotes the cut between halo and disk stars in the 13h field (which is why stars can be found "outside the boundaries": the color 
cut is different in each field, see Fig.Q, the dashed-dotted line indicates the cutoff at (6 — r)corr = 1.8 due to the metallicity 
correction. The different symbols refer to the different fields: squares - 1 h field, triangles - 9h field, dots - 13h field, open 
circles - 16 h field, crosses - 23 h field. 



Table 2. The scaleheights of the thin disk, assuming a scale- 
length of hr = 2.0kpc) and hi{hr — 3.5kpc), respectively. 
The numbers in brackets are the corresponding values of x^. 



CADIS field 


hi{hr = 2.0kpc) 


hi{hr — 3.5kpc) 


Ih 


249 ± 25 (0.92) 


255 ± 26 (0.92) 


9fi 


336 ±40 (1.21) 


312 ±35 (1.21) 


13h 


249 ± 24 (0.62) 


255 ± 25 (0.63) 


16fi 


321 ±33 (0.41) 


319 ±32 (0.41) 


23 h 


319 ± 33 (1.23) 


317 ±33 (1.21) 



hr = 2.0 kpc we find that the weighted mean of the scaleheight 
of the thin disk is hi = 281 ± 13, for hr = 3.5 kpc we find 
hi = 283 ± 13. These values are indistinguishable from each 
other, and also the values of indicate that the value adopted 
for hr does not influence the measurement of the scaleheight 
hi. 

We then fitted equation (jSj to the data over the range 
< z < 4 kpc, again for hr = 2.0 kpc, and hr — 3.5 kpc, re- 
spectively, keeping the corresponding values of hi fixed. With 
hi and (rti ±712) known, the free parameters are now the scale- 
height h2 of the thick disk, and the ratio n2/{ni ±712). Again 
the measurement is insensible to the scalelength assumed in 
equation (|5j- This is due to the high Galactic latitude of the 



CADIS fields (45 < 6 < 60) - the halo becomes dominant 
already at a projected radial distance of only about 7 kpc. 

Although the number of degrees of freedom is already 
kept to a minimum, this fit is highly degenerated, as can be 
seen from the distribution (see Fig. 13, where we asumed 
hr = 3.5 kpc. The contour lines show the 1,2 and 3cr confi- 
dence levels. Neither the ratio 712/(711 + 7*2) nor the value of 
the scaleheight of the thick disk can be estimated with high ac- 
curacy. Except for the 9 h field, where the distribution sug- 
gests a very large value of h2, the minimum can be found at 
around 900 ^ ft.2 ^ llOOkpc, and 712/(771 + 772) ^ 0.05. 
The fit of the scaleheight of the thick disk in the 9h field 
is mainly influenced by the one data point at z sa 3.6 kpc, 
where the density seems to be extraordinarily high. Because 
this value appears to be an exception, we disregard it in the cal- 
culation of the mean. The mean x^, excluding the 9h field, is 
/i2 ~ 1000, 712/(711 + 712) « 0.04. Including the 9h field, the 
mean (lower right corner) is approximately at /12 ~ 1200, 
712/(771 ±712) ~ 0.04. 

Figure Is] shows the density distribution of the disk stars in 
the five fields, and the corresponding best fits. The value for 
the local normalisation is given by the CNS4, and the values 
for hi can be found in Table |2] (we assumed hr = 3.5 kpc). 
The fraction of thick disk stars in the local normalisation is 
assumed to be 112/(711 + 712) ~ 0.04, and the mean value of 
/i2 = 1.0 kpc was adopted for all fields except for the 9 h field, 




Fig. 4. plane for the parameters /i2 and n2l(n\ + n^) for the fit of the measured density distribution in the five fields, and the 
mean of all five ones (bottom right). The scalelength is assumed to be = 3.5 kpc, and the values of are given in Table 
121 The contour lines show the \a,2a and 3a confidence levels. 
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Table 3. The exponent a from equation (|5}, for a flattened halo 
with an axis ratio of c/a = 0.6, and for c/a — 1.0 (no flatten- 
ing). The numbers in brackets are the corresponding values. 



CADIS field 


a(c/a = 0.6) 


a{c/a = 1.0) 


Ih 


2.36 ±0.14(1.37) 


2.8 ±0.16(1.05) 


9h 


0.44 ±0.08(13.99) 


0.5 ±0.1(13.68) 


13h 


1.50 ±0.09(3.13) 


1.98 ±0.09(1.54) 


16h 


2.84 ±0.14(3.24) 


3.52 ±0.18(1.65) 


23 h 


2.98 ±0.18(11.81) 


3.68 ±0.22(7.01) 



where a formal value of ft,2 = 2.6 kpc is required to achieve a 
satisfying fit. 

4.2. The stellar density distribution of the Galactic 
halo 

For the local normalisation we take stars from the CNS4, which 
have been discrimir iated against disk stars by their metallici- 
ties and kinematics ('Fu chs&.Tahreilil998!) . From these we se- 
lected stars with —0.2 < {b — r) < 0.7 (which corresponds to 
3.5 < Mv < 7.5). All these stars have very low metallicities 
(\Fe/H] < — Idex) and large space velocities (see Table 1 in 
iFuch s & Jahr eiB ( 1998)), which clearly identifies them as halo 
stars. The following six stars in a radius of r < 23 pc around the 
sun satisfy the criteria: Gl 53A, Gl 451 A, Gl 158A, GJ 1064A, 
GJ1064B, LHS 2815. Thus we find po = 89125 ± 35650 
stars/kpc^. 

Fig. ^ shows the completeness corrected density distri- 
bution of all stars. We fitted equation (|5} to the last four data 
points, while keeping the normalisation at z = fixed to 
the value given by the stars from the CNS4. We estimated 
the power law index a for a flattened halo with an axis ra- 
tio of c/a — 0.6 (dashed line in Fig. (|6j), and no flattening 
(c/a = 1.0, solid line). The values of a we deduced from the 
fit to the data are listed in Table |3l 

Clearly the fit in the 9 h field is dominated by one excep- 
tionally high data point at z w 10 kpc, and does not give a 
reliable result. The weighted mean, disregarding the 9 h field, 
is a = 2.12 ± 0.06 for (c/a) = 0.6, and a = 2.51 ± 0.07 for 
(c/a) = 1.0, where the values favour a halo without strong 
flattening: The sum of the values of the fits in the four single 
fields which enter in the determination are J2x^ — 19.55 for 
the flattened, and J^X^ — 11-25 for the spherical halo. 

The axis ratio assumed for the estimation of the power 
law index a influences the fitted value: the flatter the halo, 
the smaller a. But not only the axial ratio influences the 
measurement of a, also the local normalisation adopted for 
the fit. In order to estimate the effect we fitted the data for 
^uphm ^ fj^^ p{,°*i'"^ = Po - CTpo, respectively, 

and found weighted means of q;(pq'^''"^) = 2.26 ± 0.05 (with 
'Zx^ ^ 10.82) and a(p;)°"'"") = 1.68±0.07 (X) = 51.29) 
for (c/a) = 0.6, anda(poP'™) = 2.95 ± 0.07 (X) = 16-45) 
anda(/9[,°''"'°) = 1.97±0.07(^ = 30.12) for (c/a) = 1.0, 
respectively (again leaving out the 9 h field). 
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The values of a are getting larger for a larger value of po, 
and smaller for smaller po- However, the the summed indi- 
cate that the normalisation given by the CNS4 stars is not too 
far away from the "true" normalisation. The true value is defi- 
nitely not smaller If it is slightly larger, then a flattened halo is 
favoured. However, it is not possible to measure the flattening 
of the halo accurately with the current data set. 

4.3. The stellar luminosity function of the disk 

Knowing the density distribution of the thin disk stars, we can 
calculate the stellar luminosity function $(My). We restrict 
the analysis to stars within a radius of r < 1.5 kpc, for beyond 
this range the contribution by thick disk and halo stars becomes 
dominant. As in Paper I, we calculate absolute visual magni- 
tudes from the R magnitudes to make comparison with litera- 
ture easier. In the relevant magnitude intervall (5 < My < 14, 
i.e— 1 < (b—r) < 1.8) there holds an empirical linear relation: 

Mvj = 1.058 -Mr^. (12) 

For every luminosity bin the effective volume is calculated by 
integrating the distribution function along the line of sight, r, 
where the integration limits are given by the minimum between 
1.5 kpc and the distance modulus derived for upper and lower 
limiting apparent magnitude: 

-Rmax 

V^l^^LU J y{r,b)rUr, (13) 

-Rmin 

where 

n _ -|nO-2(16"'"'-Mi,)-2.0 
^niiii — J 

i?,„ax = min(1.5kpc,10"-2(23»--A/,)-2.0)_ 
The distribution function 

z^(r, 6) = exp(— r sin fe//ii) (14) 

is normalised to unity at z = 0. 

The weighted mean of the SLFs of the five CADIS fields is 
shown in Fig.0 For each single field we also calculated the ef- 
fective volume by using the measured values of hi. Also shown 
in Fig.0are the SLFs one would measure for hi = 0.250 kpc 
and hi = 0.350 kpc (dotted and dashed lines in the upper 
panel), in order to demonstrate how the SLF changes with the 
scaleheight assumed for the calculation of the effective volume. 
The dash-dotted line in the middle shows our old determination 
as publishe d in Paper I (Fig. 10), the dashed line is the local 
SLF (.Tah reiR& Wielenlll997l) of stars inside a radius of 20 pc 
around the sun, which is based on HIPPARCOS parallaxes. In 
the range where the SLFs overlap, they are consistent with each 
other. As in Paper I we calculated the weighted mean of CADIS 
and local SLF, which is shown in the lower panel of Fig.0 in 
comparison with a photometric SLF which is based on HST 
observations dCiould etalJl998l) . 

We will use this luminosity function as input for the simu- 
lation, which we describe in the following section. 
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5. A Monte Carlo model of the Milky Way 

The method described in the previous section has the disadvan- 
tage that it involves a completeness correction, since faint stars 
can only be traced at small distances. However, it is possible 
to avoid completeness corrections, by modeling the observed 
colors and apparent magnitudes of the stars in the fields, using 
Monte Carlo simulations. The best-fitting model of the Galaxy 
can be determined by x^-htting or a maximum likelihood anal- 
ysis. 

The two basic incredients are the stellar density distribution 
function and the stellar luminosity function. We will keep the 
luminosity function fixed, while varying the structure parame- 
ters of the density distribution. 

5. 1. The Density Distribution 

We use the same model of the Milky Way as described in 
Section Q: a sum of two exponential disks with scaleheights 
hi < h2, and a power law halo. 

The local density is again given by the CNS4 stars 

(OiireiB & Wielen 1997; Fuchs & JahreiB 1998). The contribu- 
tion of thick disk stars to the local normalisation, n2/('^i +^2), 
is left as a free parameter, and is varied from 2% to 1 8% in steps 
of 1 %. The scaleheights hi and ft,2 are varied from 150 pc to 
400 pc in steps of 50 pc, and from 200 pc to 1800 pc in steps of 
100 pc, respectively. The simulation was run for a radial scale- 
length hr = 2000 pc, and for hr = 3500 pc. The power law in- 
dex of the stellar halo was increased from a = 2.0 to a — 4.0 
in steps of 0.5, for a spherical as well as for a flattened halo 
with an axis ratio (c/a) = 0.6. 

For each field, the expected number of stars N in the model 
under consideration is calculated by integrating the distribution 
function along the line of sight: 

/•oo 

N= p{r,l,b)ujr^dr , (15) 
Jo 

where oj = 7.6145 • 10"^ for one CADIS field and p{r, I, b) = 
p{r, I, 6)disk + pir, I, &)haio- For each of the 24480 combina- 
tions of structure parameters we simulated 1000 times more 
stars than observed in each CADIS field, in order to avoid large 
poisson noise in the model data. 

5.2. Magnitudes and colors 

In order to simulate "observed" apparent magnitudes and col- 
ors and to keep the two different approaches consistent, we 
have essentially to invert the steps described in Section|3 First 
of all each star which has been simulated according to the den- 
sity distribution function is assigned an absolute magnitude 
following the stellar luminosity function determined from our 
sample (see Section l43l . The faint end of the luminosity func- 
tion of halo stars ha s only been determined with high uncer- 
tainties llGouldl2003h . Although it has been suggested that the 
luminosity function of t he halo stars might be different from 
the one of the disk stars jRevle & RobinlEoOli) . there is no ac- 
curate measurement available: The measurement of the lumi- 
nosity function of halo subdwarfs is extremely sensible to the 



distance range in which it is determined - the SLFs measured 
from local subdwarfs (where the number of stars is very small 
and limits the accuracy of the measurement) differ from those 
measured in the outer halo. The reason for these discrepancies 
might be the insufficient accuracy with whi ch the stellar struc- 
ture of the halo is known iDigbv et alJiool . Therefore, as long 
as the the luminosity function is not determined with higher ac- 
curacy and the density distribution of the halo is known better, 
we use our measured SLF for both disk and halo stars. 

In the range 6 < My < 13 the mean of the measured 
CADIS and the local SLF from|Iahrei6 & Wielen (1997) can 
be described quite precisely by a third-order polynomial in log- 
space: 

log$(Afv) = ^ qM^ , (16) 

j=0,3 

with Co = 4.78676, ci = -2.63075, C2 = 0.297266, 
-0.0104039. 

Each star is assigned an absolute V band luminosity, ran- 
domly chosen from a distribution following equation M6\ . 
Since there is no V or comparable filter available in the CADIS 
filter set, we have to convert the My magnitudes into M^^ 
magnitudes (see Paper I), i.e. Mr^, 0.945 • My. CADIS 
{h — r)c colors can be calculated from the absolute Rc band 
magnitudes using the main sequence relation described in 
Paper 1. The same relation was used in Section |4] to deduce 
absolute magnitudes from the (6 — r)c colors of the stars. 

5.2.1 . Metallicities 

In Section |4] we have corrected for the influence of different 
metallicities on the color-magnitude relation. Now we have to 
simulate those effects in order to compare the model with the 
data. Halo stars are again assumed to be fainter by AM^j = 
0.75 as disk stars of the same color In the direct determina- 
tion we have not been able to separate thin from thick disk 
stars. In the simulation however, we can account for the fact 
that thick disk stars have subsolar metallicities, and hence have 
a different color-magnitude relation than thin disk stars with 
solar metallicities. From a sample of 5 thick disk stars with 
-0.3 < [Fe/H] < -0.95, where HIPPARCOS paraUaxes, 
spectroscopic metallicities and photometry in blue and red 
passbands are available, we estimate the offset from the solar 
main sequence valid for the thin disk stars to be AMjf « 0.35 
(C. Flynn, priv.comm.). 

In order to take the spread in metallicities into account, the 
simulated colors are then smeared out by adding a value taken 
from a uniform distribution of numbers in the range —0.2 < 
(5(6 -r) < +0.2. 

5.2.2. Photometric errors 

Our observed data have of course photometric errors, which 
are of the order of = 0.01 in case of the brighter {R < 21) 
and cT/f = 0.15 for the faintest (R « 23) stars. These have 
to be simulated as well, in order to allow for a reliable com- 
parison with the observations. Using the distance modulus 
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m — M — 'o log r — 5, the apparent magnitudes of the stars 
are calculated. Each star is assigned a photometric error which 
is chosen randomly out of a gaussian error distribution with a 
variance corresponding to its apparent magnitude. The error is 
added to the apparent R magnitude. The procedure is repeated 
for the B magnitudes, which have been calculated from the 
simulated [h — r)c colors. The final colors are then redeter- 
mined from the new, modified magnitudes. 

5.3. "Observing" the Model 

Each simulated star has now an apparent Rc magnitude and a 
(h — r)c color, where metallicity effects and photometric er- 
rors are incorporated in the simulation. These data can now be 
compared with the observed data. 

The simulated stars are selected by their apparent magni- 
tudes (15.5 < Rc < 23). We bin the data, observed and sim- 
ulated, into a grid with steps of /S.Rc = 0.5 in magnitude and 
A(& - r)c = 0.2 in color. 

Due to our survey selection criteria the bins at bright appar- 
ent magnitudes contain only very few stars, which makes a 
fit of the model to the data extremely difficult. Instead of using 

as an estimate of the goodness of fit, we us e the C-stati stic 
for sparse sampling, which was developed bv lCashI d 19791) . It 
is based on the likelihood ratio, and is suited for small number 
statistics. In the limit of large numbers it converges to y^. 

If rii is the number of stars in the bin, the expectation 
value predicted by the model and N the number of bins, then 




For each of the 24 480 models of each field we calculate C 
by summing over the bins in the magnitude-color grid. 

The overall C is lowest for a ~ 2.5, (c/a) = 0.6, hr ~ 
3500 pc, and hi = 300 pc, while C is only marginally larger if 
(c/a) = 1.0 (and a — 3.0). There is no significant difference 
between models with hr — 2000 pc or hj. — 3500 pc. This 
is consistent with the results from the previous section. Figure 
(|8|i shows the C planes for /12 and ^2/ [ni + 712) (with h^ — 
3500 pc, hi = 300 pc, (c/a) = 0.6, and a fixed at 2.5), for 
each of the five individual fields and the mean of all of them. 

From the mean of the five fields we find the minimum of 
the C statistic to be at /12 = 900 pc and n^jini + n-i) « 0.09. 
While the value of /12 is consistent with the value derived from 
the direct measurement, the relative normalisation is higher 
However, the fit is highly degenerated, and the minimum in the 
distribution of C is streched along the 712/(^1 +n2)-axis. Thus 
the relative normalisation is only measurable with extremely 
large errors, and the rather high value of ri2/(ni + ri2) = 0.09 
should be regarded as an upper limit. 

Figure (|9j shows the color-coded representation of the den- 
sity of the stars in the magnitude-color diagram for the five 
fields. The images have been rebinned in steps of 0.1 in both 
color and magnitude. The contour plots overlaid are the simu- 
lated best-fit cases of the model {hr — 3500 pc, hi = 300 pc, 
h2 = 900pc,n2/(ni + n2) = 0.09, (c/a) = 0.6 and a = 2.5). 



The prominent structure at the left of each image consists 
of disk stars, whereas the halo stars are mainly distributed at 
the right. The simulation matches the data very well except for 
the 13 h field, where the gap between halo and disk is not mod- 
elled. In all cases the halo does not match the data very well, 
a possible explanation might be that the real halo luminosity 
function, which we do not know, is in fact significantly dif- 
ferent from the disk luminosity function which we use for the 
simulation. 

Figure (I10> shows the distribution of absolute V band mag- 
nitudes of the CADIS data (calculated from the [h — r)c col- 
ors as described in Section 15. 2> in comparison with the sim- 
ulated data, for the best fit cases. Both histograms have been 
normalised to the total number of stars, in order that they can 
be compared directly. 

With exception of the CADIS 23 h field the distribution of 
absolute V band magnitudes is matched very well by the sim- 
ulation. In the 9 h field the model can not account for the large 
peak at M y ~ 8. It is not impr obable that the Galactic halo 
is clumpy jNewberg et alJEooil) . so the the feature may well 
be an overdensity of stars with respect to the smooth powerlaw 
halo. 

In the 23 h field the number of bright halo stars is overpre- 
dicted. Since the 16 h and the 23 h fields are located at the same 
Galactic longitude and |6i6/i| ~ 1 623/1!, both fields should show 
the same distribution of color or absolute magnitudes, respec- 
tively. However, as both distributions differ from each other at 
the bright end, it seems more likely that the 23 h field is 'miss- 
ing' some stars, rather than that the model (which matches the 
16 h field more closely) overpredicts stars which are not sup- 
posed to be there. 

6. Summary and Discussion 

Our model of the Galaxy consists of three components: a thin 
and a thick disk (described by a sum of two exponentials), 
and a power law halo. The local normalisation of disk and 
halo stars has been fixe d by stars selected from the CNS4 
dJahreiB & Wielenlll997h by their colors and apparent magni- 
tudes. 

We have estimated the structure parameters of the stellar 
density distribution of the Milky Way using two different, com- 
plementary methods: from a direct measurement of the den- 
sity distribution of stars perpendicular to the Galactic plane 
{Baconian ansatz), and by modelling the observed colors and 
apparent magnitudes of the stars and comparing them to the 
data {Kartesian ansatz). 

Both methods have their advantages and limitations: the di- 
rect determination of the density involves completeness cor- 
rections because faint stars are only traced to small distances. 
These rely on the assumption that the stellar luminosity func- 
tion (SLF) does not change with distance and is the same for 
disk and halo stars. The same assumption of course is used in 
the simulation. However, the simulation avoids completeness 
corrections and thus the errors do not depend on an iterative 
correction. A fitting of a simulated color-magnitude distri- 
bution of stars works very well for a large number of stars; 
with the CADIS data the method is at the edge of feasibihty. 
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Fig. 8. The C planes of the five fields for the parameter /12 and 71,2/(^11 +"-2) (with hr = 3500pc, hi — 300 pc, (c/a) = 0.6, and 
a = 2.5). The plot in the lower right corner is the mean of the five indidual fields. The contour lines show the la, 2cr, and 3a 
confidence levels. 
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We therefore used the C statistic developed bv lCashl(ll979l) for 
sparse sampling to find the best fitting parameter combination. 

From the fit to the direct measurement of the density 
distribution of the stars, we find hi — 281 ± 13 pc if we as- 
sume hr = 2000 pc, and hi — 283 ± 13 pc if hr — 3500 pc. 
Obviously, due to the high Galactic latitude of the fields (45 < 
b < 60), we are not able to measure the scalelength (the halo 
becomes dominant already at a projected radial distance of only 
about 7 kpc). The fit of the scaleheight of the thick disk and the 
local percentage of thick disk stars is highly degenerate. The 
best estimate is ft,2 ~ 1000 kpc, and 122/ {nl + 712) ~ 0.04. 
We also fitted the density distribution in the halo and find 
a = 2.51 ± 0.07 for (c/a) = 1.0, and a = 2.12 ± 0.06 if 
we allow for a flattening of the stellar halo with an axial ratio 
of (c/a) = 0.6. The values of formally favour the spherical 
halo, but this result is not significant, since the power law slope 
and the axial ratio are highly correlated. Furthermore the mea- 
sured values of the exponent a, and the axial ratio c/a depend 
on the local normalisation assumed for the fit: slightly higher 
(lower) values of po yield slightly higher (lower) values of a, 
and a flattened (spherical) halo is favoured. 

From the Cartesian method we derive h2 ~ 900 kpc, 
n2/(nl + 722) « 0.09. We also find that the overall or 
C, respectively, is smallest for hi = 300 pc and an exponent 
a = 3.0. If we allow for a flattening of the halo with an axial 
ratio of (c/a) = 0.6, we find that a slightly smaller value of 
a — 2.5 fits the data better. The corresponding values of 
or C, respectively, are only marginally smaller for a flattened 
halo. The values of a are generally slightly larger than in the 
direct measurement. A possible explanation for this difference 
is that when we fit the simulated distribution of stars in the 
color-magnitude diagram, we fit all components of the galaxy 
simultaneously (whereas in the direct measurement we fit disk 
and halo separately), and there might be a degeneracy between 
halo and thick disk parameters. 

However, in general both methods yield essentially the 
same result, which we regard as a corroboration that both work 
correctly within their limitations. 

6. 1. Comparison with otiier results 

A very detailed compil ation of different re sults from different 
authors can be found in LSiegeletal.l(l2002h . The values for the 
scaleheight of the thick disk and the relative normalisation of 
thick disk stars range from 900 kpc to 2400 kpc, and from 1 % 
to 6 %, respectively. 

IChen et aP (l200lh carried out an investigation of Galactic 
structure using 5.8 • 10^ stars brighter than g' = 21™^^ from 
the Sloan Digital Sky Survey (SDSS), covering a total area of 
279 They presented models of the Milky Way for a large 
number of free parameters, that is, they determine by means of 
a maximum likelyhood analysis the scaleheights of the thin and 
thick disk, respectively, the relative normalisation of both disk 
components in the solar neighbourhood, the offset of the sun 
above the Galactic plane, the exponent a of the power law halo 
and the flattening of the halo. 



They find hi = 330 ± 3pc, 580 < /i2 < 750pc, 0.065 <, 
^^2/ ( ni + n?) ^ 0-13, an d a = 2.5 ± 0.3, and c/a = 0.55 ± 
0.06. lLemon et also found a significant flattening of 

the halo (c/a = 0.56 ±0.01). 

While our results lie well within the range found by most 
authors, our measurements for the scaleheights of the disk are 
quite dissimilar from the Sloan results. 

A possible explanation might be that their determination, 
the results of which differ significantly from all other determi- 
nations, suffers from the degeneracy due to the large number 
of degrees of freedom in their simulation: their fitting routine 
ma y have found a sec ondary minimum. 

iRobin et"al] 1I2OO O) found that there is a significant degen- 
eracy between the power law index a and the axis ratio of the 
halo. Their best fit to the data yields a — 2.44 and c/a = 0.76, 
but a flatter spheroid with c/a — 0.6 with a = 2 is not excluded 
either. We find the same degeneracy in our data: If we assume 
(c/a) = 0.6, we find a = 2.12 ± 0.06 (or a = 2.5 in the simu- 
lation), whereas for (c/a) = 1.0 we find that a = 2.51 ± 0.07 
( a = 3.0) fits the data better. 

6.2. Future prospects 

This investigation shows quite clearly that the fit of a sum of 
two exponentials is highly degenerate if the relative normalisa- 
tion n2/ni or the scaleheights hi and /i2 is not known. Also 
the measurement of a, the power law index of the stellar halo, 
and its axial ratio, (c/ a), are correlated and depend on the local 
normalisation. 

Therefore it certainly depends on the assumptions and the 
number of the degrees of freedom which values of hi, h2, ni, 
n2, a, and (c/a) are determined in different analyses. The more 
structure parameters are known and kept fixed, the narrower 
is the (or C) distribution. Thus it is important to combine 
results from different methods. 

The scaleheight of the thick disk for example will be deter- 
minable with high accuracy if the local normalisation of thick 
disk stars is known better. The local density of stars belonging 
to different components of the Galaxy has to be determined by 
means of kinematics and metalicities of the stars in order to 
break the degeneracy of the fit of a double exponential. 

The degeneracy between the power law slope and the flat- 
tening of the halo can be broken by simultaneously fitting den- 
sity distributions estimated from deep star counts in a large 
number of fields at different Galactic latitudes and longitudes, 
respectively. 

When the "size" and exact shape of the stellar halo is 
known, it will become much easier to determine the structure 
parameters of the thick disk. 

CADIS, although designed as an extragalactic survey, pro- 
vided a sufficient number of stars (and fields) to investigate the 
structure of the Milky Way. However, not only larger statistics, 
but also a much better knowledge of the luminosity function of 
halo stars is required to increase the accuracy of the results. 
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Fig. 1. The distribution of observed (fe— r)c colors in the fields. 
The dotted line shows the color cut between halo and disk. Note 
that the ordinate has not the same scale in all cases, whereas 
the dotted line extends to 20 in each plot. The maxima are not 
located at the same colors, but depend on the direction of the 
field. 
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Fig. 2. The stars in the five fields in a two color diagram ( (6 — r ) 
versus (r — i)). The spread becomes significant at (b—r) > 1.0. 
The solid line is the mean main sequence track, onto which the 
{b — r) colours of the stars are projected along the R isophotes 
of varying metallicity (dotted lines). For some extremely red 
stars, photometric errors are plotted. 
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Fig. 5. Density distribution of the stars in the disk. The solid Fig. 6. Density distribution of all stars. The solid lines are the 
line is the best fit for a sum of two exponentials. best fits for the disk (as in Figure (|5}) and a spherical halo, the 

dashed line is the fit for a flattened halo with axial ratio (c/a) = 
0.6. The corresponding exponents a are given in Table (|3})- 
The dotted line is the sum of disk and spherical halo (which is 
slightly favoured by the fit). 
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Fig. 7. Upper panel: the mean SLF of the five CADIS fields 
(solid line). For the calculation of the effective volumes we 
used the fitted values of h,\ (see Table |3. The dotted and 
dashed lines are the determination for h\ = 0.250kpc, and 
for h\ — 0.350 kpc, respectively. Middle panel: The solid line 
is our current measurement, as in the upper panel. The dotted 
line is our first measur ement as pu blished in Paper 1, the dashed 
line is the local SLF ( JahreiB & Wielen,.1997.) . which is based 
on HIPPARCOS parallaxes. Lower panel: weighted mean of 
CADIS and local SLF, in compariso n with a ph otometric SLF 
which is based on HST observations ( iGould et''al..l998.) . 





Fig. 10. Comparison of absolute V band magnitudes of the sim- 
ulation (dotted line) with the CADIS data (solid line), for the 
best fit case. 



